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Cold inelastic collisions of atoms or molecules are analyzed using very general arguments. In free 
space, the deactivation rate can be enhanced or suppressed together with the scattering length of the 
corresponding elastic collision via a Feshbach resonance, and by interference of deactivation of the 
closed and open channels. In reduced dimensional geometries, the deactivation rate decreases with 
decreasing collision energy and does not increase with resonant elastic scattering length. This has 
broad implications; e.g., stabilization of molecules in a strongly confining two-dimensional optical 
lattice, since collisional decay of the highly vibrationally excited states due to inelastic collisions is 
suppressed. The relation of our results with those based on the Lieb-Liniger model are addressed. 
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Feshbach resonances P, have been used to control 
atomic interactions in trapped ultracold quantum gases 
by tuning a magnetic field near a diatomic molecule Fes- 
hbach resonance to convert atoms into weakly bound 
molecules. For fermionic atoms the molecules formed 
were remarkably long-lived [3] , whereas for bosonic atoms 
in a BEC , collisional decay of the highly excited vibra- 
tional molecular state occurs and only a small fraction 
of molecules is observed in this case. 

Here we show, using general scattering theory argu- 
ments, that inelastic ultra-cold collisions in reduced di- 
mension can be strongly suppressed. This generalizes the 
work of Ref. which uses the exactly solvable Lieb- 
Liniger many-body model for indistinguishable bosons in 
one-dimension (ID) to all kinds of quasi-lD scatter- 
ing processes occurring in atomic waveguides, e.g., colli- 
sions of atoms and molecules. Quasi-ID scattering occurs 
in a gas in the presence of a waveguide potential that con- 
fines a 3D gas sufficiently tightly in two directions so the 
radial confinement energy luj_ (in units where h = 1) is 
much larger than the collision energy Id , as in 2D optical 
lattices |9j, elongated atomic traps '10^, and atomic in- 
tegrated optics devices 11). This suppression has broad 
implications, e.g., it can be used to stabilize molecules 
produced from bosonic atoms in tight atomic waveg- 
uides, since vibrational-to-translational energy-transfer 
collision rates are significantly reduced relative to the 
3D rates at low collision energy. Suppression of inelas- 
tic scattering can also occur in collisions of other excited 
collision partners (e.g., in hyperfine excited atom colli- 
sions) . Long-lived excited ultracold atoms and molecules 
may be useful in quantum interferometry and quantum 
computation. 

The theoretical framework for calculating atom-diatom 
scattering can be drawn along the lines of the Arthurs 



and Dalgarno model The scattering state j^*) can 

be expressed in term of a sum over basis functions. 



|*)^^^,(r)|x,). 
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where r is the atom-diatom relative coordinate, V'i(i") is 
the relative wave function, and \xj) includes molecular 
and center-of-mass degrees of freedom for channel j. We 
shall not require details of this model since our argu- 
ments are very general, and in fact extend beyond this 
particular problem (e.g., molecule-molecule collisions). 

Low-energy inelastic exoergic collisions in the presence 
of a Feshbach resonance can often be treated as multi- 
channel scattering with zero-range interactions described 
by boundary conditions for s-wave radial wave functions 
^jir) = ^/ ipj{r)di}r (see [l^ and references therein). 
In our case of low energy inelastic scattering resulting 
in deactivation of the excited state of a molecule, the 
boundary conditions take the form. 
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for the input channel (po, the closed channel (pc, and the 
deactivation product tpd having a set of output chan- 
nels {d}. This method is applicable to collisions of any 
type of particles when s-wave scattering is allowed. Note 
that broad Feshbach molecules j3| cannot be treated us- 
ing the zero-range approach of Eq. ||2Jl. However, if one 
is considering atom-molecule or molecule-molecule colli- 
sions, the resonance does not coincide with the resonance 
in atom-atom collisions, e.g., two resonances in collisions 
of Cs9 molecules have been observed at 12.72 and 13.15 
G [13, far off the atom-atom resonance at 19.84 G. 

When the coupling of the input channel to the other 
channels vanishes, Eq. reduces to the Bethe-Peierls 
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boundary condition, 



dr 



r=0 



flbg 



-ipo{0) {uncoupled} , (3) 



where Obg is the non-resonant background elastic scatter- 
ing length. Hence, from 0, we find that C/qo = — l/obg- 
Outside the interaction region, V'o(r) satisfies the 
Schrodinger equation. 



2M 



V2^o(r) + yconf(r)Vo(r) = EM^^), (4) 



where l^onf is the confining waveguide trapping potential, 
E is the collision energy, and M is the reduced mass of the 
colliding particles. Moreover, the radial wave functions 
ifc and ipd satisfy the Schrodinger equations 



1 d'^ipcM 

2M dr^ 
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where Dc is the asymptotic value of the closed channel 
potential, is the deactivation energy for channel d, 
and we assumed that Vconf ^ Dc.d- These equations can 
be solved to obtain 



= fc{0)cxp i^-^2M {D,-E): 
(fidir) = Rdexp{ipdr) , 



(6) 
(7) 



where pd = \/2M [E + Dd). The closed channel has an 
attractive potential (C/cc < 0) and a single bound state 
with energy E^^^^ ^ - U^J (2M). 

Substitution of Eqs. ® and ^ into 101 leads to the 
following boundary condition: 



dipo{r) 



dr 



¥^0(0). 

flcff 
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The deactivation energies typically substantially exceed 
all interaction energies. Therefore only the contributions 
of zero and first orders in \Ujk\/pd need be retained. To 
this accuracy we can neglect terms proportional to Udd' 
in |(5J), so the amplitudes Rd can be expressed as 
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with 



abg|t/< 
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Here, A = cn,g\Ucc\\Uoc\'^ / {fJ-M) is the resonance 
strength, fi is the difference of the magnetic moments 
in the closed and open channels, and 



Es - 



\Ucc\ 

M 



(9) 



is an effective detuning. 



The length UcS has an imaginary part due to coupling 
to the deactivation channels and can be expressed as 



Ocff — Obg 



Es - iVc 
Es+iiA-iT' 



where 



Pd 



abg C/oc 



Uo 



(10) 



— ahglUdol'^Es 



For a tightly bound closed-channel state, or when the 
detuning of the collision energy from the Feshbach energy 
is small, lE'posh - < U^^/M, then Es « i^Fosh - E. In 
this case, neglecting deactivation, one can approximate 
flcff by the real effective energy-dependent length , 



Ocff {E) K, flbg 



1 
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E^^i{B-Bo)\ 
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where B — Bq = A + -Eposh//^ is the detuning of the 
external magnetic field B from its resonant value Bq. 
The deactivation cross section can be expressed as 



d 



(12) 



where ipo corresponds to the input channel wave function 
normalized to unit incident flux density, and the factor 



S 



^Pd 



Udo 
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(13) 



describes interference of deactivation of the closed and 
open channel states. 

First, we consider collisions in free space {Vconi = 0). 
The proper solution of Eq. Q has the form 



V Po 



exp (ipo • r) — 
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tpo r 



- exp {ipor) 



(14) 

where the collision momentum is po = V2ME. For low- 
energy collisions, when abgPo ^ 1, we find 



M 
V Po 

The deactivation cross-section 

AttS. 



Cfr, 



PO 
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(15) 



(16) 



diverges at low collision energies, while the deactivation 
rate coefficient 



AttS. 



(17) 



has a finite non-zero limit proportional to loefip. The 
deactivation is suppressed at i?^ = 0, where Oeff is close to 
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zero. Under certain conditions it can also be suppressed 
due to interference in the factor S. The deactivation is 
enhanced near resonance of Ceff at Es = — /iA. In this 
case it can reach a maximum value of 



rGe,inax 



'W 



(18) 



Consider now collisions in a transverse harmonic 
waveguide. This problem has been analyzed in Refs. |^ 
Ir^ for a single-channel Huang pseudopotential, which is 
equivalent to the boundary condition in Eq. The 
case of a multichannel (5-function interaction has been 
considered in Ref. using a renormalization procedure. 
Equations (17) and (19) in ^ express the proper solu- 
tion of Eq. Q in terms of the transverse Hamiltonian 
eigenf unctions |nO) with zero angular momentum projec- 
tion on the waveguide axis, z, 



V'o(r) 




exp(ipo2;) 1 00) 



-Ma 



conf 



exp {ipn\z\) \nO) 
n=Q \Jn- (poa±/2)^ 



(19) 



Here a± = (Mw_l) ' is the transverse harmonic oscil- 
lator length, u>±_ is the waveguide transverse frequency, 
Pn - ^2M[E-{2n+l)uji_], 



Tconf (Po) 



Ma I 



(20) 



is the transition matrix, and C,{v, a) is the Hurwitz zeta 
function 0, . The wave function H19|) is normalized 
so the average incident flux density per waveguide area 
7ra^ is unity. The sum in Eq. H19|l diverges as r ^ 0. 
The divergent part can be evaluated as a^/r 0|. This 
leads to V3o(0) = - \ Ma^ i/M/poT'conf {po ) , and to the 
deactivation rate coefficient 



Krnuf = vrAfa , IT, 



conf I 



(21) 



For weak confinement, when aj^pa 3> 1, approximation 
(49) in Ref. 0| leads again to Eq. (fnji for the wave 
function and to Eq. H17(l for the deactivation rate. For 
strong confinement, when a±pQ <C 1, approximation (41) 
in Ref. [l^ leads to 



Tconf(Po) ~ -i 



M 



-Cai_pQ - I- 



(22) 



where C ~ 1.4603. At low collision energies, or at large 
Ocff, where po ^ |aoff|/alj the wave function at the ori- 



gin, v?o(0) 



lj_^/Mpo, is much less than the corre- 



sponding value of fo{0) in free space given in Eq. (fT5|l . 
Thus confinement prevents the particles from occupy- 
ing the same position. A similar effect is responsible for 




FIG. 1: (color online) Scaled deactivation rate coefficient 
as a function of scaled magnetic field detuning b = 



fi{B — Bo)/{2uj±) — i and collision energy e = E/{2uj±) — i, 
calculated for Obg = 0.1a± and = 2uj±. In the confined ge- 
ometry, the solid and dashed curves correspond to e = 10~^ 
and e = 0.5 respectively, whereas the free space results are 
given by the dot-dashed and dotted curves respectively. The 
inset shows deactivation rate versus e in confined geometry 
for 6 = (solid curve) and b = 100 (dashed curve) . 



fermionization of ID bosons with strong interactions 
Under these conditions the deactivation rate. 
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can be substantially suppressed by confinement. 

This conclusion is graphically demonstrated in Fig. ^ 
under conditions when Oeff is expressed by Eq. Hll|) . It 
shows resonances in the deactivation rate at Es — —fiA 
for collision energies comparable to u}± and in free space, 
as well as deactivation suppression near Es — 0. At low 
collision energies, when 



Po < kbgl/a? 



(24) 



deactivation under confinement does not have resonances 
and can be strongly suppressed even compared to the 
non-resonant process in free space. Suppression appears 
also at £' = (2fc + l)uj±, where excitations of transverse 
waveguide modes become o pen, leading to jumps in the 
elastic scattering amplitude pi3.ll3|. 

The above results are obtained for a system composed 
of two arbitrary types of particles interacting via s-wave 
scattering. A suppression of inelastic collision has been 
predicted in Ref. |y] for a many-body system of ID in- 
distinguishable bosons using the Lieb-Liniger model j7[. 
However, as we shall see below, the suppression is mostly 
a two-body interaction effect even in this model. 

Consider first the two-body scattering process with 
particle momenta pi and p2- The Lieb-Liniger wave func- 
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tion with unit norm in interval [0, L] (i — > oo) is, 



, P1-P2 



1 



V2L 

imUa 



Pi -P2 



imU„ 



exp (ipizi + ip2Z2 
exp (ipiZ2 + ip2Zi) 



(2) 

for < zi < Z2 < i and ^'p/pa (-^1,22) 



(2) 



(^2,^l) 



for < Z2 < 



- P1P2 V-'i I -^2; — ^piP2 
zi < L. Here m is the particle mass, 

and Ua ~ 4abg [ma^(l — Cabg/a^)] is the interaction 
strength ■sj . The two-body correlation function with the 
particles at the same position, 



92^ {Pl,P2) 



(0,0) 



{Pi -P2) 



L2 



(25) 

is the probability to find two particles at the same place. 
Equation (|25|l already describes qualitatively the behav- 
ior of g2 when the ratio of the interaction to collision 
energies is large, as obtained in 32 ~ {pi — P2Y /U"^- 
In the A^-body case, the two-body correlation function 
52^' can be estimated as a sum of 9^2^ over all pairs of 
the colliding particles with the quasimomenta pj and pj' , 



92\Pj^Pj') 

j<j' 



(2) 

dpidp2f{pi)f{p2)92 {Pi,P2), (26) 



(27) 



in full agreement with the results of Ref. P|. In the 
Tonks-Girardeau regime, f{p) ~ l/(27r) for \p\ < irp and 
f{p) = otherwise, and (|26|l leads to 



27rV 
372 



(28) 



This value is half the exact value determined in Ref. [6j . 
The difference can be related to higher order correlations, 
which are neglected here (the two-body picture can only 
include second order correlations). However, (|28|l de- 
scribes the correct behavior of §2^^ as 7 — > 00, leading 
to suppression of all kinds of collision phenomena under 



tight confinement when mUa/p >> 1 [this condition has 
the same meaning as Eq. (|24ll ]. 

In summary, inelastic collision rates in free space 
demonstrate resonances and dips, being proportional to 
Iflefip, and are capped by H18|l . Interference of deactiva- 
tion of the open and closed channels can also suppress the 
rate. In quasi-lD scattering at low collision energies [see 
Eq. (|24() ]. inelastic collisions do not have resonances and 
are suppressed. This result applies to the collision of any 
type of atoms or molecules interacting via s-waves and is 
not based upon the ID Bose gas Lieb-Liniger model. 
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